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^ ■ Abstract 

O . 

The f(R) gravity models formulated in Einstein conformal frame are equivalent to 
Q ' Einstein gravity together with a minimally coupled scalar field. The scalar field couples 

■ with the matter sector and the coupling term is given by the conformal factor. We 

use this interacting model to derive a necessary condition for alleviating the coincidence 
Q ' problem. 
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There are strong observational evidences that the expansion of the universe is accelerating (see 



e.g. [1]). However, the origin of this cosmic acceleration is not well understood and remains 
as one of the main challenges of modern cosmology. The standard explanation invokes an 
O ■ unknown component, usually referred to as dark energy. It contributes to energy density of 
£sj ■ the universe with = 0.7 where fid is the corresponding density parameter [2]. A candidate 
for dark energy which seems to be both natural and consistent with observations is the cosmo- 
logical constant [2] [3] [4] . However, in order to avoid theoretical problems [3] , other scenarios 
have been investigated. In one of these scenarios the matter sector remains unchanged and the 
gravitational part suffers from some modifications. A family of these modified gravity models 
is obtained by replacing the Ricci scalar R in the usual Einstein-Hilbert Lagrangian density 
for some function f(R) [5]. 

There are two important problems that are related to the cosmological constant. The first 
problem, usually known as the fine tuning problem, is the large discrepancy between obser- 
vations and theoretical predictions on its value. There have been many attempts trying to 
resolve this problem [3]. Most of them are based on the belief that the cosmological constant 
may not have such an extremely small value at all times and there should exist a dynami- 
cal mechanism working during evolution of the universe which provides a cancelation of the 
vacuum energy density at late times [6]. The second problem concerns with the coincidence 
between the observed vacuum energy density and the current matter density. While these 
two energy components evolve differently as the universe expands, their contributions to total 
energy density of the universe in the present epoch are the same order of magnitude. Besides 
the possibility that the present epoch may be a stationary regime at which the ratio of the two 
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energy densities are constant, it is also quite possible that we live in a very special epoch, a 
transient epoch at which the ratio varies slowly with respect to the expansion of the universe. 
A possible solution to the coincidence problem is to consider an interaction between dark en- 
ergy and dark matter. If such an interaction exists the two corresponding energy densities do 
not scale independently. It is shown that, this can lead to a constant ratio of energy densities 
when an appropriate coupling term is applied [7] [8]. 

In the present note, we will consider the coincidence problem in Einstein frame representation 
of f(R) gravity models. In these models the dynamical variable of the vacuum sector is the 
metric tensor and the corresponding field equations are fourth order. This dynamical variable 
can be replaced by a new pair which consists of a conformally rescaled metric and a scalar 
partner. Moreover, in terms of the new set of variables the field equations are those of Gen- 
eral Relativity. The original set of variables is commonly called Jordan conformal frame and 
the transformed set whose dynamics is described by Einstein field equations is called Einstein 
conformal frame. The dynamical equivalence of Jordan and Einstein conformal frames does 
not generally imply that they are also physically equivalent. In fact it is shown that some 
physical systems can be differently interpreted in different conformal frames [9] [10]. The 
physical status of the two conformal frames is an open question which we are not going to 
address here. Our motivation to work in Einstein conformal frame is that in this frame there 
is a coupling between the scalar degree of freedom and matter sector induced by the confor- 
mal transformation. As previously stated, there is a large amount of interest to realize the 
coincidence problem as a consequence of an interaction between matter systems and the dark 
sector. Although the whole idea seems to be promising, however, the suggested interaction 
terms are usually phenomenological and are not generated by a fundamental theory. In our 
case the interaction term is given by the conformal factor. We investigate the consequences of 
this interaction term and derive an expression which constrains the form of the f(R) function. 
We will show that this constraint selects those f(R) models that allow for possible alleviation 
of the coincidence problem. 

2 Framework 

The action for an f(R) gravity theory in the Jordan frame is given by 



where k = 8ttG, G is the gravitational constant, g is the determinant of g^ v and S m is the 
action of (dark) matter which depends on the metric and some (dark) matter field ip. 
Stability in matter sector (the Dolgov-Kawasaki instability [11]) imposes some conditions on 
the functional form of f(R) models. These conditions require that the first and the second 
derivatives of f(R) function with respect to the Ricci scalar R should be positive definite. The 
positivity of the first derivative ensures that the scalar degree of freedom is not tachyonic and 
positivity of the second derivative tells us that graviton is not a ghost. 

It is well-known that f(R) models are equivalent to a scalar field minimally coupled to gravity 
with an appropriate potential function. In fact, we may use a new set of variables 
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where Q = = /'(-R) and /3 = This is indeed a conformal transformation which 

transforms the above action in the Jordan frame to the following action in the Einstein frame 
[9] [12] 

Sef = \S d ' x ^ - V ^ - 2 W)} + S m (9^e 2 ^, V) (4) 

All indices are raised and lowered by g^ v unless stated otherwise. In the Einstein frame, (f> is 
a minimally coupled scalar field with a self-interacting potential which is given by 

Vim) = ZkMR) (5) 

Note that the conformal transformation induces the coupling of the scalar field with the 
matter sector. The strength of this coupling j3, is fixed to be and is the same for all 
types of matter fields. In the action (4), we take g^ v and as two independent field variables 
and variations of the action yield the corresponding dynamical field equations. Variation with 
respect to the metric tensor g^, leads to 

GV = fc(7)t + Tp (6) 

where 

1%, = V^V„0 - i^V^0V 7 - V{<t>)g, v (7) 

are stress-tensors of the scalar field and the matter field system. The trace of (6) is 

V 7 0V 7 + 4V(<j>) - R/k = f m (9) 

which differentially relates the trace of the matter stress-tensor T m = g^T™ to R. Variation 
of the action (4) with respect to the scalar field 0, gives 

□0-^M = -/3v^T™ (10) 
d(p 

It is important to note that the two stress-tensors T™ v and Tfi v are not separately conserved. 
Instead they satisfy the following equations 

V"7^ = -V"T^ = f3VkV v (p T m (11) 

We apply the field equations in a spatially flat homogeneous and isotropic cosmology described 
by Friedmann-Robertson- Walker spacetime 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) (12) 
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where a(t) is the scale factor. To do this, we take T™ and T£ v as the stress-tensors of a 
pressureless perfect fluid with energy density p m , and a perfect fluid with energy density 
P<p = \4> 2 + V{4>) and pressure = \4> 2 — V ((/>), respectively. In this case, (6) and (10) take 
the form ' 

3H 2 = k( P<t> + p m ) (13) 
2H + 3H 2 = -kuj^ (14) 

+ 3#0 + = -(3Vkp m (15) 

d(p 

where = ^ is equation of state parameter of the scalar field 0, and overdot indicates 
differentiation with respect to cosmic time t. The trace equation (9) and the conservation 
equations (11) give, respectively, 

2 + R/k - 4V(<f>) = p m (16) 

p m + 3Hp m = Q (17) 

P(f) + 3H(cu <p + 1) P4> = -Q (18) 

where 

Q = f3Vk<j>p m (19) 

is the interaction term. This term vanishes only for = const., which due to (3) happens 
when f(R) linearly depends on R. The direction of energy transfer depends on the sign of Q 
or 0. For > 0, the energy transfer is from dark energy to dark matter and for < the 
reverse is true. 

We emphasize that the coupling term (19) is very similar to some phenomenological coupling 
terms suggested in the literature. In fact, there are different kinds of interacting models which 
have been investigated [7] [8] . A particular class of these models considers Q = acpp in which a 
is a coupling constant, ip is usually a quintessence field and p is energy density of dark matter 
[8]. Apart from the similarity of the latter with (19), there are also some important differences. 
Firstly, the scalar field is not a kind of matter field and is actually given in terms of the func- 
tion f(R). Secondly, j3 is a universal coupling constant implying that couples with the same 
strength to all types of matter fields. In contrary, it is possible to consider a as a non-universal 
coupling constant so that it may couple to dark matter and baryons with different strengths 
[13]. Moreover, the value of (5 is fixed to be l/v^6 while a is constrained by observations [8]. 
We will return to this last point later. 



3 The Coincidence Problem 

One of the important features of the cosmological constant problem is the present coincidence 
between dark energy and dark matter energy densities [14]. There is a class of models in which 
this observation is related to some kinds of interaction between the two components [7] [8]. In 
these models the two components are not separately conserved and there is a flow of energy 

^Hereafter we will use unbarred characters in the Einstein frame. 
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from dark energy to dark matter or vice versa. In this sense, dark energy and dark matter 
energy densities may have the same scaling at late times due to the interaction, although 
they decrease with the expansion of the universe at different rates. The important task in 
this context is to find a constant ratio of dark energy to dark matter energy densities for 
an appropriate interaction term. Despite the fact that this approach seems to be promising, 
there is not still a compelling form of interaction which is introduced by a fundamental theory. 
Therefore one usually uses different interaction terms and tries to adapt them with recent 
observations. 

In f(R) gravity models presented in the Einstein frame, there is a fixed interaction between 
the scalar field and matter sector. Since the form of the interaction is fixed by the conformal 
transformation one can therefore search for some appropriate forms of the function f(R) for 
which the energy densities ratio of the two components takes a stationary value. This is the 
strategy that we are going to pursue in this section, namely, to find some conditions on the 
functional form of f(R) that may lead to a constant r = p m / Pcj, ■ 
To do this, we consider time evolution of the ratio r , 

Pm P<f> /on\ 

r = r— (20) 

P<i> P<p 

From equations (17), (18) and (19) we obtain 

f = 3#rc^ + @Vk<j)r(r + 1) (21) 

In this relation, we can write r in terms of the parameters r and q. We first use (13) and (14) 
to replace the equation of state parameter oj<f> with the deceleration parameter q. Applying 

H = -(q + l)H 2 (22) 

to the equation (14) gives 

(2g - 1)# 2 

UJ<p - 7= (A5J 

v kp^ 

We then use (13) in the latter and substitute the result in (21), which leads to 

f = Hr(2q - l)(r + 1) + pVkjrr(r + 1) (24) 

On the other hand, we can combine the trace equation (16) with the equations (5) and (13) 
to obtain 

^ = i { ?rr + 3i/2(2,; - 3)(1 -7 ) - 2 ?j } (25) 

When we put this expression into (24), the result is an equation that relates r to the parameters 
r, q and H. The requirement that the universe approach a stationary stage in which r either 
becomes a constant or varies more slowly than the scale factor, leads to the following relation 

g(f';H,r s ,q) = (26) 
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where 

3r 2 2f i 

tf, r s , g) = r s (2q - l)(r s + 1) + /3r s (r s + Vi^fj + 3(2g - 3)(1 - -) - (27) 

and r s is the value of r when it takes a stationary value. It is now possible to use (26) to 
check that whether a particular f(R) model is consistent with a late-time stationary ratio of 
energy densities. In general, to find such f(R) gravity models one may start with a particular 
f{R) function in the action (1) and solve the corresponding field equations for finding the 
form of q(z) or H(z). However, this approach is not efficient in view of complexity of the field 
equations. An alternative approach is to start from the best fit parametrization q(z) obtained 
directly from data and use this q(z) for a particular f(R) function in (26). Here we will follow 
the latter approach. 

For a given redshift Zq and the parameters r s (zo), q(zo) and H(zq), the relation (26) acts as 
a constraint on the function f(R). As an illustration, we apply this constraint to some f(R) 
functions. Before doing this, there are some remarks to do with respect to (26). This condition 
is a consequence of r = when r = r s becomes stationary at late-times. At sufficiently late- 
times characterized by z = zo, we take r s = r and rewrite (26) as 

g(f o \;H o ,r Q ,q o ) = (28) 

where 

^;^ ,ro,go)=ro(2go-l)(ro + l)+/3ro(ro + l){^ T + 3(2g -3)(l-^)-2-^}3 (29) 

r o + 1 Jo -"o/o 

Here the functions / , f' Q and f ' are the late-time configurations of f(R), f'(R) and f"(R) 
which are obtained by replacing R with 

R = Q(l-q)H 2 (30) 

at the redshift zq. Note that an f(R) gravity model is usually given in terms of some pa- 
rameterizations. In this sense, the condition (26) acts actually as a constraint relating the 
corresponding parameters of a particular f(R) gravity model to the constants q , r and H . 
We use a two-parametric reconstruction function for characterizing q(z) [15] [16], 

1 012 + 92 , Q1 x 

q{z) = 2 + JT+zf (31) 

Fitting this model to the Gold data set gives q ± = 1.47±i;ll and q 2 = -1.46 ± 0.43 [16]. We 
also take z = 0.25 which, with use of (31), corresponds to q ps —0.2. Moreover, recent 
observations imply that r = P ^f z ^ ~ f [17]. 
Now let us first consider the model [18] [19] 

f(R) = R + XRoA n (32) 

Here R is taken to be of the order of Hq and A, n are constant parameters. In terms of the 
values attributed to these parameters, the model (32) is divided by three cases [19]. Firstly, 
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when n > 1 there is a stable matter-dominated era which does not follow by an asymptoti- 
cally accelerated regime. In this case, n = 2 corresponds to Starobinsky's inflation and the 
accelerated phase exists in the asymptotic past rather than in the future. Secondly, when 
< n < 1 there is a stable matter-dominated era followed by an accelerated phase only for 
A < 0. Finally, in the case that n < there is no accelerated and matter-dominated phases 
for A > and A < 0, respectively. Thus the model (32) is cosmologically viable in the regions 
of the parameters space which is given by A < and < n < 1. 

When we use (30) in the function g^fa H ,r ,q ), it takes the form of an expression which 
relates the parameters n and A to go, r o and H . In fig.l we have plotted g(n, A; H , r , go) for 
A = —1. This figure indicates that the constraint (28) is satisfied only for n « 0.9 which im- 




0.70 0.75 0.80 0.85 0.90 0.95 

n 



1.00 



Figure 1: The plot of g(n, A; H Q ,r ,q ) for the model (32) when A 
r = 3/7. The vertical dashed line corresponds to n = 0.906. 
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-0.2 and 



plies that for this value of the parameter n, the model (32) admits a late-time stationary ratio 
of the energy densities. Note that n f» 0.9 lies in the range that the model is cosmologically 
viable. 

Now we consider the model presented by Starobinsky [20] [21] 



f(R)=R- 1 R {l 



[l + (^) 2 ]" m } 



(33) 



where 7, m are positive constants and Ro is again of the order of the presently observed 
effective cosmological constant. Using the same procedure, we have plotted the function 
g(m, 7; Ho, ro, go) in fig. 2. The figure shows that there are some regions in the parameters 
space for which the condition (28) is satisfied. The condition is satisfied on the upper bound- 
ary of the surface plotted in fig. 2 where g(m, 7; H , r , go) = 0. Thus for the corresponding 
values of the parameters, the coincidence problem can be addressed in the context of the model 
(33). For instance, as the figure indicates the parameters space is bounded by 7 > 10.5 and 
m > 0.04 so that for m > 0.04 the parameter 7 should remain near the value 10.5. 
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Figure 2: The plot of g(m, 7; H , r , go) for the model (33) when q = —0.2 and r = 3/7. 

4 Conclusion 

In Einstein frame representation of f(R) gravity models, the scalar partner of the metric 
tensor interacts with (dark) matter in such a way that the interaction term is fixed by the 
conformal transformation. This means that contributions of the scalar field and the (dark) 
matter system to total energy density do not scale independently. As a consequence, even 
tough the two components may start with different scalings at early times, they may have the 
same scaling at sufficiently late times. 

We have considered this feature as a possibility for addressing the coincidence problem. In 
fact, the interaction of dark energy and dark matter has been recently taken as a natural 
guidance for alleviating the coincidence problem by some authors. In absence of an interaction 
or coupling term based on a fundamental theory, most of the current investigations have been 
limited to a phenomenological level. In our case, the interaction term, Q, is given by the 
conformal transformation and can be written in terms of R, f'(R) and f"(R). Due to stability 
considerations, any viable f(R) model should satisfy f'(R) > and f"(R) > [22]. Thus 
the direction of the energy transfer is determined by the sign of R in a particular epoch. For 
instance, in an epoch for which R > 0, the energy transfer is from dark energy (or the scalar 
field 4>) to dark matter while for R < the reverse is true. 

We have derived a relation giving the evolution of the parameter r. We have found that there 
is a class of f(R) gravity models satisfying the condition (28) for which a late-time stationary 
state for r exists. As illustrations, we have shown that the model (32) lies in this class only 
for n ~ 0.9. The condition is also used for the Starobinsky's model. We have shown that 
there is a region in the parameters space for which the coincidence problem can be addressed 
in this model. The region is characterized by the upper border of the surface plot of the fig. 2 
for which g(m, 7; H , r , q ) = 0. 

Finally, we point out that there is no a free parameter in the interaction term (19) since /3 
is fixed by conformal transformation. In general, the interaction of the scalar field and the 
matter sector may lead to a fifth force and violation of equivalence principle. In fact, the real 
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challenge for alleviating the coincidence problem comes from the combination of restrictions 
from local gravity experiments and dynamical considerations. Thus the question is that how 
a coupling term without a free parameter can be consistent with local gravity experiments. 
The point is that, in our case, these experiments constrain the corresponding parameters of a 
particular f(R) gravity model 1 '" rather than the coupling constant of the interaction term. For 
the model (33), it is shown [24] that the most stringent bound is m > 0.9 which comes from 
violation of equivalence principle. Combining the latter with the bounds indicated in fig. 2, one 
infers that alleviation of the coincidence problem requires that 7 ~ 10.5. 
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